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Abstract 



A program package, which facilitates computations in the framework of Analytic 
approach to QCD, is developed and described in details. The package includes the 
explicit expressions for relevant spectral functions calculated up to the four-loop 
level and the subroutines for necessary integrals. 
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1 Introduction 



The strong interactions display two fundamental features, namely, the asymptotic 
freedom (at high energies) and color confinement (at low energies). The first 
feature allows one to study the strong interaction processes in the ultraviolet 
domain by making use of perturbation theory. However, theoretical description 
of hadron dynamics at low energies requires nonperturbative methods. 

In general, there is a variety of nonperturbative approaches to handle the 
strong interaction processes at low energies. In this work we will focus on 
the so-called "dispersive" (or "analytic" ) approach to Quantum Chromodynam- 
ics (QCD). The basic idea of this approach is to merge the perturbative results 
with the nonperturbative constraints arising from relevant dispersion relations. 
In turn, this eliminates some intrinsic difficulties of perturbation theory and ex- 
tends its range of applicability towards the infrared domain. One of the pos- 
sible implementations of this approach within QCD is the so-called "massless" 
Analytic Perturbation Theory [TH3]. The latter has been successfully employed 
in the studies of various strong interaction processes (see, e.g., papers [U-tH], 
reviews [21 [151 [16] an d references therein). The incorporation of effects due to 
the nonvanishing mass of lightest hadron state has been implemented into ap- 
proach in hand within the so-called "massive" Analytic Perturbation Theory, see 
Refs. [HUTS] for the details. 

A central object of the current approach is the so-called spectral function, 
which can be calculated by making use of the strong running coupling a s (Q 2 ) 
(see Sect. [3]). At the one-loop level the perturbative running coupling and rele- 
vant spectral function have a quite simple form (Eqs. (132|) and (Q, respectively). 
However, at the higher loop levels the strong running coupling has a rather cum- 
bersome structure, and the calculation of corresponding spectral functional rep- 
resents a rather complicated task. 

The primary objective of this paper is to calculate (by hands) the explicit 
expressions for the aforementioned spectral functions at the higher loop levels and 
to incorporate them (together with proper subroutines for necessary integrals) 
into a single program package. In turn, the latter will facilitate computations 
within the approach in hand. 

The layout of the paper is as follows. Section [2] constitutes a brief description 
of the Analytic approach to QCD. The calculation of spectral functions and the 
numerical evaluation of basic integrals is considered in Sects. [3] and 0J respec- 
tively. Section describes the template program and includes the list of package 

: At the higher loop levels the spectral functions can also be calculated numerically. However, 
it requires a lot of computation resources and essentially slows down the overall computation 
process. 
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commands. The basic results are summarized in the Conclusions (Sect. [6]). Ap- 
pendix |A] contains explicit expressions for the perturbative strong running cou- 
pling up to the four-loop level. Appendix |B] contains explicit expressions for the 
spectral functions calculated up to the four-loop level. 



2 Analytic approach to QCD 

In general, in the framework of perturbation theory the high energy behavior of 
the strong correction d(Q 2 ) to a physical observable D(Q 2 ) can be approximated 
by the power series in the strong running coupling a s (Q 2 ). Namely, at the £-loop 
level 

dfUQ 2 ) = E^K } (Q 2 )] J = S>(^Y \^\Q 2 )t <? -+ °°' « 

3=1 i=i 

where Q 2 = —q 2 > is the spacelike kinematic variable, dj stands for the rele- 
vant perturbative expansion coefficient, a^ e \Q 2 ) is the £-loop perturbative strong 
running coupling (see App. [SJ, f3 = 11 — 2n f /3 denotes the one-loop perturba- 
tive (3 function expansion coefficient, ri{ is the number of active quarks, and 
a^(Q 2 ) = a^(Q 2 )/3o/(4vr) is the so-called "couplant". However, perturbative 
expansion ([T]) is valid in the ultraviolet domain only. In particular, at any loop 
level dpe Tt (Q 2 ) (OQ) possesses unphysical singularities in the infrared domain. This 
fact contradicts the general principles of local Quantum Field Theory and essen- 
tially complicates the analysis of low energy experimental data. Besides, pertur- 
bative expansion ((T|) can not be directly employed in the theoretical description 
of physical observables depending on the timelikd^] kinematic variable s = q 2 > 0. 

As it has been noted in the Introduction, one can overcome the aforementioned 
difficulties of perturbative approach by invoking relevant dispersion relations. 
Thus, in the framework of "massless" Analytic Perturbation Theory (APT) [THS] 
the theoretical expressions for the strong corrections d(Q 2 ) and r(s) to physical 
observables D{Q 2 ) and R(s), depending on spacelike (Q 2 = —q 2 > 0) and timelike 
(s = q 2 > 0) kinematic variables, take the form (see papers [2J[2] and references 
therein for the details) 

d^(Q 2 ) = td J {^jA^ J ( Z ), Afi J (z) = ] g f^±da, z = %, (2) 

3 

rS T ( S )=E^(^j4tiH> A^( W ) = J g f\a) d ^ v ) = ±. (3) 

J— 1 " w 



2 For example, the experimental data on the so-called i?(s)-ratio of electron-positron an- 
nihilation into hadrons can only be examined by making use of both perturbation theory and 
relevant dispersion relation, see, e.g., Ref. [T§] . 
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Here and further A denotes the QCD scale parameter and Qj (er) stands for the 
spectral function corresponding to j-th power of the £-loop couplant (see Sect. E]). 

In general, the effects due to the mass of the lightest hadron state can be safely 
neglected at intermediate and high energies only. In particular, such effects play 
an essential role in theoretical description of the strong interaction processes at 
low energies. As it has been noted in the Introduction, the effects due to the 
nonvanishing mass of the lightest hadron state have been accounted for within 
so-called "massive" Analytic Perturbation Theory (MAPT) [T71IT8] . Thus, in the 
framework of MAPT the theoretical expressions for the above-mentioned strong 
corrections read 

d^(<?,m^=f i d i (^) j AU j (z, x ), z = %, X = ^, (4) 

4?,(,,)^/> M ^, (5) 

x 

r$ APT (s,m 2 ) = J2 d j{^) 42,jO,x), w = -j^, (6) 
A^ J (w, X ) = 9(w- X )J f\a)^, (7) 

w 

where 9(x) is the unit step function (9(x) — 1 if x > and 9(x) =0 otherwise). 
It is worth noting that in the limit m — ?> Eqs. (TjO) and ([6]) coincide with Eqs. (ED 
and ([3]), respectively: g?m1p T (<5 2 , 0) = d!il T (Q 2 ) and rMl PT (s,0) = rip T (s). Be- 
sides, 4!iip T (Q 2 ,m 2 ) -»■ dfl T (Q 2 ) for Q 2 > m 2 and r^l PT (s,m 2 ) = r£. T (s) for 
s > m 2 (see Refs. [T7T.[T8] for the details). It is worthwhile to emphasize also that 
expressions (j2J), (JSJ) and (SI) — (JTJ) represent the nonpower functional expansions^ of 
the strong corrections d{Q 2 ) and r(s). 



3 Calculation of the spectral functions 



In general, there is no unique way to restore the aforementioned spectral func- 
tion g~j (a) by making use of the perturbative expression for the strong running 



coupling a^(Q 2 ) (discussion of this issue can be found in Refs. 
what follows we adopt the definition proposed in Refs. [TH3]: 

1 



). In 



lim 

2ni e-K)+ 




-A 2 



IE) 



(0 



-A> 



IE) 



A 



(f) 



SL,l( Z ) 



3 In other words, for j > 2 Az) ^ 

less, in the ultraviolet asymptotic (z — > oo) A$l j(z) 
and (|3]) reproduce the perturbative power series ([lj. 



andA^-H^ A%[ 1 (w) 



TL, j 



Nonethe- 



{Q 2 )\ 3 ', i.e., the expansions ([2]) 
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Figure 1: The one-loop perturbative couplant (Eq. fl32|) . dot-dashed curve) and 
the one-loop first-order expansion functions: APT (Eq. (jTUl) . dashed curves) and 
MAPT (Eqs. ([II]), (H2D, solid curves). The values of parameters: A = 350 MeV, 
m = 270 MeV. 



(cr is a dimensionless variable). At the one-loop level (£ = 1) the first-order (j = 1) 
spectral function ([H]) can easily be calculated (see Eq. 



— lim fa« 

2ni e->o+ 



lim 

2ni £->o+ 



-A 2 (a + ie] 



-A 2 (<x - zel 



1 



?/ 2 + 7T 2 



(9) 



ln(— a — ie) \n(—a + ie) 

where y = In a. Eventually, this leads to the following expressions^ for the one 
loop (£ — 1) first-order (j = 1) expansion functions (j2J), ©, (JSJ), and (J7|): 



i i 

+ 



In 2! ' 1 — 2 

1 2 1 + X 



4l,iW 



1 1 / In w " 
arctan 

2 7T V 7T 



In z + x + <x + z 



a 



1 1 / In w " 
arctan 

2 7T V 7T 



(10) 

(11) 
(12) 



see papers [SHE] and references therein for the details. The functions (1T0~ ]) - (IT21 
are shown in Fig. [TJ 



4 It is assumed that arctanx is a continuously increasing function of its argument: —ir/2 < 
arctan a; < n/2 for — oo < x < oo. 
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Figure 2: The four-loop spectral functions Qj (cr) (jHi): f?i (c) (solid curve), 
£>2 ( <T ) (dotted curve), 10 ■ 03 4 (cr) (dashed curve), and 10 2 • q£*((t) (dot-dashed 
curve). The values of parameters: rif — 3, y = In a. 



At the higher loop levels (£ > 2) perturbative couplants a^(Q 2 ) have a cum- 
bersome structure (see App. [A}. Hence, the calculation of the £-loop spectral 
function of j-th order gf (cr) ((8]) (1 < j < £) represents a rather complicated 
task. The numerical calculation of the spectral functions (jSJ) requires a lot of 
computational resources and essentially slows down the running of the program. 
Nonetheless, this problem can be resolvec^ in the following way. 

First of all, it is convenient to express the spectral functions gfp (a) (jSJ) in terms 
of the real and imaginary parts of the £-loop couplant a^(Q 2 ) on a physical cut: 



lim 

£->0+ S 



-A 2 (a =f ie) 



= A 



(«), 



(13) 



In this case Eq. (jSJ) can be represented in a concise form: 



y( 3 )r 1 )k 2k 
t \2k + V 



2fc+l 



42W 



j-2k-l 



(14) 



where Kj = j/2 + (j mod 2)/2 — 1, and 



n 



,my ml (n — m)\ 



(15) 



5 An alternative way to overcome this problem is to construct a set of explicit expressions 
which approximate the nonpower expansion functions ([2]), ([3]), ([5]), and (|7J) accurately enough, 
see, e.g., Ref. [2T] . 
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is the binomial coefficient. In particular, at any loop level 



9f {a) 



4» 



m 



(16) 
(17) 



— 7T 



7T 



(19) 



Then, one has to calculate (by hands) the real and imaginary parts ( |T3l) of 
the couplants a^(Q 2 ) (j32)-(|35) on a physical cut. In turn, this will enable 



one to construct the spectral functions gf\cr) 



14]) corresponding to any integer 
power (j > 1) of the couplant a^(Q 2 ) up to the four-loop level. The explicit 



expressions for the calculated functions A^(a) and A^(a) (1 < £ < 4) are given 
in App. [B] (see also App. C of Ref. [6] and App. C of Ref. [16]). The four-loop 
[I = 4) spectral functions Qj{&) (1 < J ' < 4) are shown in Fig. [21 
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4 Basic integrals 



The explicit expressions for the spectral functions gf\cr) (fl4"|) enable one to 
compute the corresponding nonpower expansion functions in spacelike (Eqs. (121) 
and ([5])) and timelike (Eqs. ([3]) and (J7])) domains. Specifically, in the framework 
of massless APT the spacelike expansion functions ([2]) read 



AsJz) 



a + z 



da. 



(20) 



For the numerical evaluation of A SL (z) it is convenient to split Eq. (12"U|) into three 
terma^l, namely 



r(x)e l l x 
e i/x + ~ 



dx + 



Qyiv) 
1 + ze~y 



dy + 



r(x) 



1 + ze- 1 ^ 



dx. 



(21) 



where 



Qyiv) = 



a=ev 



Qx{x) = Q y {y) 



y=l/x 



r(x) 



(22) 



3 This was first proposed by Dr. I.L. Solovtsov. 
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Figure 3: The four-loop "massless" APT expansion functions A s 2 j{z) © and 
^rl j{ w ) © : j = 1 (solid curves), j = 2 (dotted curves, scaled x 10), j = 3 
(dashed curves, scaled xlO 2 ), and j = 4 (dot-dashed curves, scaled xlO 3 ). The 
values of parameters: rif = 3, A = 350 MeV, z = Q 2 /A 2 , w = s/A 2 . 



The numerical integration of Eq. ( l2~Tj) is implemented within QCDMAPT library 
by the subroutine APTSL (see Sect. [5]). In turn, the timelike massless APT ex- 
pansion functions © take the form 



A Th (w) 



da 

q{<t) —. 
a 



(23) 



Similarly to the previous case, it is worth representing this equation in the fol- 
lowing way: 

1 



■h -1, 



A Th (w) = < 



hiw 

J 2 (lnio, 1) + MO, 1), 
*(<> 1 



where 



In w 



■f2 



+ J 2 (-l,l) + Ji(0,l), \nw<-l 
-1 < law < 1 

1 < lnu? 



(24) 



W2 



Ji(xi,x 2 ) = / r(x)dx 



(25) 



The numerical integration of Eq. (1241) is implemented within QCDMAPT library 
by the subroutine APTTL (see Sect. |5]). The four-loop massless APT expansion 
functions A^j(z) and A® tj (w) (1 < j < 4) are shown in Fig. [3j 
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Figure 4: The four-loop MAPT expansion functions Ag L j -(2;,x) fl5]) and 
^■TL i^jX) : J = 1 (solid curves), j = 2 (dotted curves, scaled x 10), 
j = 3 (dashed curves, scaled xlO 2 ), and j = 4 (dot-dashed curves, scaled xlO 3 ). 
The values of parameters: n f = 3, A = 350 MeV, m = 270 MeV, z = Q 2 / a2 > 
u; = s/A 2 . 



In the framework of massive APT the spacelike expansion functions ([5]) read 



A s d z , x) 



z + X 



a -X da 

Q{°) — ; • 

a + z a 



(26) 



For the numerical evaluation of A SL (z, x) it is convenient to represent Eq. (12"6"|) in 
the following way: 



z 



A SL (z,x) = < 



z + X 

z 



/l(-l,j^-)+/ a (-l,l)+/3(0,l) 



lnx < -1 



z + X 

z 



J 2 (lnx,l) + / 3 (0,1) 
1 



-1 < lnx < 1 



I z + x V ' In y/ ' 



h 0, 



1 < lnx 



where 



X2 1 /r 

x 2 ) = / r(x) — dx. 



j-i 
m 

h (2/1,2/2) = / 



e l/x + ^ 



e y - X 
ey + z 



dy, 



(27) 



(28) 
(29) 



•'- • -1/x 



l:d-r\..C->) = I / >'-'-) ~ + *l- 1/x dr ' (;J,()) 



X'l 



The numerical integration of Eq. ( )27|) is implemented within QCDMAPT library 
by the subroutine MAPTSL (see Sect. [5]). In turn, the timelike MAPT expansion 
functions (j7]) take the form 

°° da 

A TI >, x) = 0(w -x) q{°) — = 6(w- X )A tl (w), (31) 

J a 

w 

where A Th (w) is defined in Eq. (T24"|) . The numerical integration of Eq. (l3Tj) is 
implemented within QCDMAPT library by the subroutine MAPTTL (see Sect. [5]). 
The four-loop MAPT expansion functions A^j(z, x) and A^Aw, x) (1 < j < 4) 
are shown in Fig. HI 



5 Description of the template program 

The provided template program (template .mw) is self-explaining. Besides, both 
template program and QCDMAPT library (QCDMAPT. lib) contain detailed com- 
ments. 

The template program is organized as follows (see template.pdf). On the 
first page the QCDMAPT library is loaded and a sample set of input parameters 
is specified: n f = 3 active quarks, Q = 700 MeV, s = (700 Me V) 2 , m = 270 MeV, 
and A = 350 MeV. The evaluation of expansion functions for the given set of pa- 
rameters is presented on page 2. In particular, the first computation deals with 
the first-order (j = 1) £-loop (£ = 1, 2, 3, 4) expansion functions, namely, a^(Q 2 ) 

(Eqs. m-m), (Eq- ©), 4?,iH ( E( l- ©), 4£i(*,X) (Eq- ©), and 

A^l i(w, x) (Eq. ©). The second computation evaluates four-loop (£ = 4) expan- 
sion functions of j-th order (j = 1,2,3,4), namely, [a^{Q 2 )] j (Eq. (|55))), 4?,j(^) 

(Eq. ©), 4ti(^) (EQ- ©), 4?,i(^X) (Eq- ©), aud^Kx) (Eq. ©)• The 
data arrays of the first-order one-loop expansion functions are computed and rel- 
evant plots are presented on page 3 (see Fig. [1] and its caption for the details). 
The data arrays of the massless APT four-loop [I = 4) expansion functions of 
j-th order (j = 1, 2, 3, 4) are computed and corresponding plots are presented on 
page 4 (see Fig. [3] and its caption for the details). The data arrays of the MAPT 
four-loop (£ = 4) expansion functions of j-th order (j = 1, 2, 3, 4) are computed 
and relevant plots are presented on page 5 (see Fig. H] and its caption for the 
details). 

The summary of QCDMAPT package commands is given below. 
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• PERTURBATIVE QCD 

AlphaPertlL(z) = a^(Q 2 ) {£ = 1,2,3,4; z = Q 2 /A 2 ): £-\oop perturbative 
strong running coupling, see App. lAl 

APertlL(z) = af\Q 2 ) = of )(Q 2 )/3 / (4tt) {£ = 1,2,3,4; z = Q 2 /A 2 ): £-loop 
perturbative QCD couplant, see Eqs. (I3"2"|) - (j3"5|) . 

beta j — /3j (j — 0, 1, 2, 3): perturbative (j + l)-loop expansion coefficient of the 
QCD function, see Eqs. (PHP]) . 

Bj= Bj = (3j/f3^ +1 (j = 1,2,3): combination of QCD (3 function perturbative 
expansion coefficients, see App. lAl 



• SPECTRAL FUNCTIONS 



ARe£LY(y) 



4>w 



a=ev 



[l = 1,2,3,4): real part of the £-loop perturbative 



QCD couplant on a physical cut, see Eqs. f fl3|) . (122]) . and App. 



AIm£LY(y) = A^(a 



a=ev 



1,2,3,4): imaginary part of the £-loop pertur- 



bative QCD couplant on a physical cut, see Eqs. (1T3"]) . (12"2"]) . and App. 



Rho£LjpY(y) 



a=ev 



(£ = 1,2,3,4; j = 1, 



spectral function (in 



terms of y = In a variable) corresponding to j-th power of the 
turbative QCD couplant, see Eqs. (151). AH]). (1221. and App. EE 



-loop per- 



R£LjpX(x) 



rf\x) 



(£ = 1,2,3,4; j = 1, ...,£): spectral 
function (in terms of x = 1/ In a variable) corresponding to j-th power of 
the £-loop perturbative QCD couplant, see Eqs. (IE]), ( TT4|) . ( )22|) . and App. [El 



• BASIC INTEGRALS 



APTSL(RhoY,RX,z) = A SL (z): numerical integration for the massless APT space- 
like expansion functions, see Eqs. (120]) -( 122]) with g y {y) =RhoY(y) and 
r(x) =RX(x). 

APTTL(RhoY,RX,w) = A Tlj (w): numerical integration for the massless APT time- 
like expansion functions, see Eqs. (|23|) - (|25|) . (|22l) with g y (y) =RhoY(y) and 
r(x) =RX(x). 

MAPTSL (RhoY , RX , z , chi ) = A st (z, x)- numerical integration for the MAPT space- 
like expansion functions, see Eqs. (I26l)-(l30l). ( 1221) with g y (y) =RhoY(y) and 
r(x) =RX(x). 
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MAPTTL(RhoY,RX,w,chi) = A TL (w,x)- numerical integration for the MAPT 
timelike expansion functions, see Eqs. (l3Tj) . (T23|) — (T25|) . fl22|) with g y (y) = 
RhoY(y) and r(x) =RX(x). 

•APT EXPANSION FUNCTIONS 

ASLlLjp(z) = A® tj (z) (£ = 1,2,3,4; j = 1, ...,£): APT spacelike expansion 
function (^-loop level, j—th. order), see Eq. (J2J). 

ATLlLjp(w) = A^l^w) {£ = 1,2,3,4; j = !,...,£): APT timelike expansion 
function (£-loop level, j—th. order), see Eq. (jSJ. 

AlphaSLlL(z) = A^ x {z) 4tt/(3 (£ = 1,2,3,4): APT spacelike £-loop effective 
coupling, see Eq. (j2J). 

AlphaTLlL(w) = A® tl (w) 4tt//3 (£ = 1,2,3,4): APT timelike Moop effective 
coupling, see Eq. (j3J). 

•MAPT EXPANSION FUNCTIONS 

ASLmlLjp(z,chi) = A$j(z,x) (t = 1,2,3,4; j = !,...,£): MAPT spacelike 
expansion function (£-loop level, j-th order), see Eq. (JSJ). 

ATLmlLjp(w,chi) = A^^w^) = 1,2,3,4; j = 1, ...,£): MAPT timelike 
expansion function (£-loop level, j-th order), see Eq. (J7J). 

AlphaSLml L(z,chi) = A$ tl (z, x) 4vr//3 (£ = 1, 2, 3, 4): MAPT spacelike £-loop 
effective coupling, see Eq. (j5J). 

AlphaTLmlL(w,chi) = A® A (w, x) 4tt//3 {£ = 1, 2, 3, 4): MAPT timelike £-loop 
effective coupling, see Eq. ([7]). 

6 Conclusions 

A program package, which facilitates computations in the framework of the An- 
alytic approach to QCD, is developed. It includes the explicit expressions for 
relevant spectral functions calculated up to the four-loop level and the subrou- 
tines for necessary integrals. 

This work was partially supported by the grants RFBR-08-01-00686, BRFBR- 
JINR-F08D-001, NS-1027.2008.2, and JINR grant. 
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A Perturbative QCD running coupling 



This Section contains explicit expressions for the perturbative strong running 
coupling af>(Q 2 ) = Att a)p (Q 2 ) / 0o up to the four-loop level. Specifically, 



af\Q 2 ) 

4 3) (Q 2 ) 



af\Q 2 ) 



1 



In z 
1 

In z 
1 

In z 



-B 1 



B 1 



Q 2 

ln(lnz) 
~kFz~' 

ln(lnz) 
In 2 z 



ln J z 



In z 



-B 1 



B\ [in 2 (hi z) -ln(ln^) - 1 
ln(lnz) 



Br 



ln 2 z 



ln 3 z 



\n 4 z 



B{ ln^(lnz) -ln(ln^) - 1 



+ B, 



B\ 



ln 3 (lnz) + -ln 2 (lnz) 
2 



(32) 
(33) 



(34) 



+21n(lnz) 



W,B 2 ln(lnz)+-5 3 



(35) 



where Bj = (3j/ is the combination of QCD j3 function perturbative expansion 
coefficients: 



0o 
0i 

02 
03 



11 n f , 

3 ' 

38 

102 - — n { , 

2857 5033 325 2 

rif H n r , 

18 54 f ' 



2 

149753 



(36) 
(37) 
(38) 



6 



3564C(3) - 



1078361 6508 



+ 



50065 6472 
L 162 + 81 ^ ' 



n\ + 



162 27 
1093 q 



C(3) 



rif 



729 



-n f . 



(39) 



In these equations rif stands for the number of active quarks and ((x) denotes the 
Riemann ( function, ((3) ~ 1.202. The one- and two-loop coefficients (/3q and /3i) 
are scheme-independent, whereas the expressions given for (3 2 and (3 3 are calcu- 
lated in the MS subtraction scheme (see papers [22TT25] and references therein for 
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the details). Note that expressions fl33l) - fl35|) correspond to approximate solutions 
of the perturbative renormalization group (RG) equation for the QCD invariant 
charge, see, e.g., Ref. [26], Sect. 2 of review [T5], and App. A of review [IB] . 
Nonetheless, the difference between the approximate running coupling a^ e \Q 2 ) 
(!33l) - (l35|) and the exact solution of perturbative RG equation a^(Q 2 ) is not 
controllable at every considered loop level. 



B Spectral functions 



As it has been mentioned in Sect. El the spectral functions Qj (<j) 
represented in the following form: 



»5V> 



£ 

k=0 



2k 



(-l) k ir 2k 



where 



£->0+ 



lim [-A 2 (a =p is)] = A$ e (a) =F 



2fc+l 



42M 



j-2k-l 



can be 



(40) 



(41) 



£-loop perturbative couplants a^(Q 2 ) are given by Eqs. fl32|) - fl35|) . and Kj = 
j/2 + (j mod 2)/2 — 1. The functions Af^a) and A^ e (a), calculated up to the 
four-loop level (1 < £ < 4), are listed below. 

One-loop level: 



y 2 + 7r 2 ' 



y 2 + 7r 2 ' 



In a. 



(42) 



Two-loop level: 



where 



(?/ 2 + 7r 2 ) 2 
1 

(y 2 + n 2 ) 2 



2ySf\y) + ^ 2 -y 2 )Sf\y) 
(y 2 -n 2 )S{ 2 \y) + 2n 2 ySi 2 \y) 



S[ 2 \y) 



2 i ^ 
7T , 



1-fii 



1 1 



arctan 



7T 



(?) 



(43) 
(44) 

(45) 
(46) 



and Bj = (3j/[3^ +1 is the combination of QCD /3 function perturbative expansion 
coefficients (see App. |A|) . 
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Three-loop level: 



(y 2 + 7T 2 ) 3 
1 



(y 2 + 7T 2 ) 3 



(3y 2 -7r 2 )4 3) (y)-y(y 2 -37r 2 )^ 3) (y) 
y(y 2 - 3vr 2 )^ 3) (y) + vr 2 (3y 2 - n 2 )S?\y) 



4%) 

4 3 eV) 

where 

Sf 3) (y) 

Sf(y) 

Sffo) = y 2 -7r 2 -B 1 [y5[%)-7r 2 ^(y) 



±ln(y 2 + 7r 2 ), 

1 1 /sA 

arctan — , 

2 7T W 



+Bf{s[%) [S[%) - l] - vr 2 [^(y)] 2 - l} + B 2 , 



Si 3 \y) = 2y-B 1 \sf\y)+ySf\y) +BlSf\y)2Sf\y)-l 



c(3) 



(3), 



5 2o(3), 



.--(3), 



Four-loop level: 
42(a) = 

4 4 e } (-) = 



4y(y 2 -7r 2 ) S?\y) + S?>(y) + Sf>(y) 



( 4 ), 



(2/ 2 + vr 2 ) 4 

+ [4vrV - (y 2 - vr 2 ) 2 ][^ 4) (y) + $%) + ^ 4) (y) 
(^t{[^ 2 - ^ " 4 -V][si 4) (y) + 4 4) (y) + Si 
+4n 2 y(y 2 - tt 2 ) [^ 4) (y) + (y) + # (</)]}, 



(4) 



where 



^ 4) (y) 
4 4) (1/) 



ln(y 2 



7T 



1 1 /sA 
5 2 y{^ 4) (y)[^ 4) (y)-l 

-5?7^>(y)[2S< 4 >(y)-l 

5i 2 fe 4) (y)^ 4) (i/)-i 



— 7T 



+ ^-1 

5? 



— 7T 



+ y^ 4) (y)[25r(y)-lj+|f-n 



7 (4) 



Bo 
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= b!\s{%) 



S?\y) 



+ 2 



(Sf"(!/)) 2 - 7T 2 (Si'\y)f 



5 9 



1 / -B? 



B 



-3 



+ 5^ 4) (y)-3f| + 2}, 

y(y 2 - 3vr 2 ) + B 1 [2ir 2 yS?\y) - (y 2 - n 2 )S ( *\y) 
3y 2 -n 2 -B 1 \( y 2 -n 2 )S?\y) + 2yS?\y) 



(59) 



(60) 

(61) 
(62) 
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